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摘要 : 射电 天 文 观测 数据 受 无 线 电 的 影响 日 益 严 重 。 因 此 ， 射 频 干 扰 ( Radio Frequency 
Interference, RFI) 消除 已 经 成 为 信号 处 理 流 程 中 不 可 或 缺 的 一 步 。 考 虑 了 自 适 应 阵列 信号 处 
理 中 的 移动 干扰 源 的 消除 问题 ， 对 于 阵列 天 线 望 远 镜 而 言 ， 射 频 干 扰 消除 可 以 通过 采样 协 方 
差 矩 阵 在 空域 实施 处 理 。 在 很 多 应 用 场景 中 ， 可 得 到 的 零点 深度 受 限 于 协 方差 矩阵 的 估计 误 
差 ， 进 而 影响 干扰 子 空 间 的 估计 精度 。 方 法 应 用 了 一 种 多 项 式 模型 以 跟踪 阵列 协 方差 矩阵 随 
时 间 的 变化 ， 消 除 干扰 子 空间 的 估计 误差 ， 提 高 干扰 消除 性 能 。 最 后 通过 仿真 对 比 了 传统 的 
子 空间 投影 (Subspace Projection, SP) 算法 和 基于 多 项 式 模 型 的 子 空间 投影 (Polynomial- 
augmented Subspace Projection, PSP) 算法, 仿真 结果 证 明了 所 提 方 法 的 有 效 性 。 
关键 词 干扰 消除 ; 自 适应 波束 形成 ; 子 空 间 投影 ; 射电 天 文 
中 图 分 类 号 : TN973 ”文献 标识 码 : A ”文章 编号 :; 1672-7673(2017)03-0297-07 


射频 干扰 是 射电 天 文中 一 个 十 分 重要 的 问题 。 大 部 分 的 射电 源 信号 属于 微弱 信号 ， 需 要 长 时 间 的 
积分 才能 够 被 接收 机 检测 。 由 于 低 质 量 的 滤波 器 泄露 和 邻近 频 融 的 互 调 和 干扰， 射电 天 文 观 测 数 据 总 受 
到 来 自 被 国际 电信 联盟 (International Telecommunication Union, ITU) 保护 的 频带 内 的 无 线 电信 号 的 干 
扰 。 男 一 方面 ， 根 据 可 观测 灵敏 度 与 带宽 之 间 的 关系 可 知 ， 最 小 灵敏 度 与 带宽 成 反比 ， 即 采用 更 宽 的 
输入 带宽 能 够 获得 更 微弱 信号 的 探测 能 力 。 那 么 和 目前 日 益 增 长 的 无 线 电 频率 使 用 必然 形成 了 矛盾， 人 
工 产生 的 比如 卫星 传输 信号 、 雷 达 信和 号、 无线 通信 信号 和 电视 广播 信号 等 大 大 强 于 射电 天 文 信号 ( 太 
阳 射电 爆发 信号 除外 ) ， 如 果 在 观测 频段 内 ， 会 造成 接收 机 饱和 、 假 谱 、 交 调 干扰 等 一 系列 问题 ， 影 
响 天 文 观测 的 数据 质量 。 

近 几 年 来 ， 一 系列 自 适应 波束 形成 算法 已 经 成 功 应 用 于 空间 干扰 消除 "了 ,但 是 仍然 有 很 多 应 用 
场景 和 信号 环境 是 目前 的 方法 不 能 应 对 的 。 当 干扰 移动 相对 较 快 时 ,传统 的 自 适 应 干扰 消除 方法 性 能 
下 降 。 本 文 方法 基于 已 被 广泛 应 用 的 子 空 间 投影 算法 ” ， 主 要 考虑 到 它 良 好 的 零 陷 深 度 。 子 空间 投 
影 算 法 是 通过 阵列 协 方差 矩阵 的 特征 值 分 解 得 到 干扰 子 空 间 的 估计 ， 继 而 将 接收 信和 号 投影 到 干扰 子 空 
间 的 正 交 空间 ， 以 实现 干扰 消除 。 

虽然 本 文 所 提 方 法 和 原理 通常 适用 于 大 范围 的 传 感 需 阵列 和 自 适 应 阵列 的 干扰 消除 ， 但 在 这 里 着 
重 强调 射电 天 文中 的 相 控 阵 天 线 (Phased Array Feeds, PAFs)。 这 个 应 用 对 干扰 移动 相对 于 阵列 协 方 
差 矩 阵 的 估计 时 间 间 隔 而 言 是 相对 快速 的 问题 ， 是 一 个 恰当 的 处 理 方法 。 采 样 时 间 间 隔 通 常 被 称 为 短 
时 积分 窗 ( Short-term integration, STI)。 

相 控 阵 天 线 提供 了 一 种 有 前 途 的 新 方法 以 利用 干扰 信号 的 空间 结构 跟踪 和 消除 射频 干扰 而 不 会 丢 
失 任何 有 用 数据 。 对 于 空间 干扰 消除 而 言 ， 上 述 方法 甚至 可 以 消除 完全 遮盖 有 用 信和 号 的 频谱 的 干 
扰 源 的 影响 。 

本 文 第 1 节 给 出 了 信号 模型 的 数学 表示 ; 第 2 节 介 绍 了 子 空间 投影 的 波束 形成 的 详细 原理 和 由 于 
干扰 源 的 相对 移动 导致 的 估计 误差 是 如 何 限制 零 陷 深度 的 ; 第 3 节 提 出 了 一 种 适用 于 移动 干扰 源 的 基 
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于 矩阵 多 项 式 模型 的 时 变 阵 列 协 方差 矩阵 ， 以 改善 传统 的 子 空 间 投影 算法 ; 第 4 节 给 出 了 仿真 性 能 分 
析 ， 由 结果 可 知 本 文 方法 的 干扰 消除 性 能 优 于 传统 的 子 空间 投影 方法 。 


1 信号 模型 


假设 x[n] 是 Px1 阶 的 阵列 天 线 的 复 采 样 序列 ， 它 可 以 表示 为 
x[n] =a,s[n] +a,[n]d[n] + nLn], (1) 
其 中 ，s[z] 为 观测 的 射电 天 文 信号 ; a, 为 s[ nj] 阵列 响应 矢量 ; d[ nj] 为 干扰 源 信 号 ; a,[n] A din] 
列 响应 矢量 。 天 线 是 固定 跟踪 观测 源 以 保持 射电 源 在 视 轴线 上 ， 因 此 a, 是 常量 。 随 着 干扰 源 的 相对 
运动 或 者 天 线 跟踪 射电 源 ，a,[n] 随 时 间 改 变 。 品 声 矢 量 m[n] 假 设 为 空域 白 高 斯 咖 声 。 信 号 s[n]、d 
[n] 和 m[n]j 假 设 为 广义 宽 平 稳 的 高 斯 随机 过 程 。 那么 ， 阵 列 协 方差 矩阵 表示 为 
R,[n] =E\x[n]x"[n]} =o°a,a" + a,l n]a} ln] + R, =R, +R,[n] + R,, (2) 
其 中 ，R,[n] 的 时 间 依 赖 性 来 自由 移动 干扰 源 产生 的 阵列 响应 矢量 的 变化 ， 因 此 在 一 个 短 时 积分 窗 
内 ， 它 只 能 近似 看 做 是 广义 平稳 的 。 这 个 假设 在 干扰 源 移动 速率 相对 于 短 时 积分 窗 宽 度 较 快 时 会 失 
bm 效 。 其 他 的 所 有 信号 假设 为 广义 平稳 的 。 
> BA R [n] Bin 连续 变化 ， 它 只 有 假设 在 第 个 短 时 积分 间隔 内 是 平稳 的 ， 因 此 可 以 得 到 它 的 工 
长 度 采样 的 估计 值 ， 


q De 


1 
R,, =— > x[n] x" [n] = —X,X7, (3) 
; L n=kL L 


Hoh, X,= [æ [kL], =, x[(k+1)L-1]] 为 第 个 短 时 积分 窗 工 长 度 的 采样 数据 。 
2 传统 的 干扰 消除 算法 


2.1 子 空间 投影 干扰 消除 

在 满足 窄带 的 条 件 下 ， 波 束 形成 器 的 输出 信号 可 以 表示 为 

yin] = w"x[n] , (4) 

RF, w 为 波束 形成 器 的 权 矢 量 ， 以 提供 期 望 信号 来 向 上 的 响应 方向 图 。 如 果 干 扰 信 号 相对 于 信和 号 在 
阵列 间 的 传播 时 间 是 时 域 宽 带 的 ， 那么 x[n] 可 以 被 分 解 为 多 个 子 带 ,每 一 个 有 各 自 的 波束 形成 权 矢量 。 
€ 在 传统 的 自 适应 波束 形成 算法 中 有 很 多 可 以 用 来 设计 w 以 得 到 期 望 的 方向 图 主因、 控制 方向 图 
© FEAR ACE TBS We) LE 。 

考虑 到 子 空间 投影 算法 可 以 得 到 更 深 的 零 陷 深 度 ， 间 投 影 算法 实现 “” | F 
扰 信 号 的 强度 大 到 可 以 主导 R,[nj] 中 的 咖 声 和 有 用 信号 时 ,波束 形成 权 矢量 可 以 用 下 面 的 方法 计算 : 


Wsp = Piw 

P, =I- U, U, 
其 中 , kh AGI ERS) a RSE Pr 为 估计 的 干扰 子 空间 的 正 交 投影 矩阵 ; U BE TIAE 
解 中 干扰 信号 对 应 的 单位 长 度 归 一 化 的 主 特征 向 量 : 


R [Ual Ung gl = Way! Ung NA. (6) 


MERE w 用 来 建立 无 干扰 情况 下 的 静态 天 线 方向 图 。 投 影 算 子 Prw 可 能 扰乱 或 者 损坏 期 望 的 
静态 方向 图 ， 因 此 需要 保证 干扰 从 天 线 方向 图 的 劳 欠 进入 。 男 一 方面 ，UV ,的 秩 必须 小 于 天 线 个 数 P 
或 有 充分 的 自由 度 形 成 远离 空域 零点 的 期 望 的 静态 方向 图 。 

在 移动 干扰 源 的 应 用 场景 中 ，ws, ,在 每 一 个 新 的 工 长 度 的 采样 中 被 重新 计算 以 保持 干扰 消除 零点 
在 干扰 来 向 上 。 在 多 干扰 源 的 应 用 场景 中 ， 必 须 事 先 计算 干扰 个 数 mw， 而 此 时 的 干扰 信号 矩阵 ,的 
阶 数 为 Pxm。 本 文 方法 着 重 应 用 于 单个 干扰 的 情况 ， 此 时 U, ,是 一 个 列 和 撩 量 。 


(5) 
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2.2 影响 干扰 消除 性 能 的 因素 

提高 采样 数 可 以 减 小 协 方差 矩阵 R, ,的 采样 估计 误差 ， 那 么 干扰 子 空间 的 估计 精度 和 干扰 消除 算 
法 的 性 能 会 随 之 提高 。 但 当 传 统 的 子 空间 投影 的 波束 形成 抗 干 扰 方 法 应 用 于 移动 干扰 源 的 场景 时 ， 如 
果 干 扰 移动 速度 过 快 而 采样 间隔 过 长 ， 会 引起 阵列 相应 的 不 一 致 ， 造 成 阵列 协 方差 矩阵 的 秩 不 等 于 干 
扰 源 数量 即 子 空 间 扩散 ， 那 么 自 适应 波束 形成 抗 干扰 算法 将 无 法 使 用 。 因 此 分 析 采 样 估计 误差 和 子 空间 
扩散 对 传统 自 适应 波束 形成 干扰 消除 算法 的 影响 ， 为 改进 传统 的 自 适应 波束 形成 干扰 消除 算法 做 准备 。 

(1) 采样 估计 误差 ， 在 平稳 信号 环境 中 太 , 是 R.[n]1,.4 的 一 致 信 计 ， 估 计 误差 的 方差 碱 少 了 了 。 
为 了 减少 估计 误差 ， 短 时 积分 窗 的 长 度 应 该 尽 可 能 长 

(2) 子 空间 扩散 : R ,包含 了 第 大 个 采样 时 间 间 隔 的 工 长 采样 的 时 变 R,[n] 的 时 间 平 均 。 尽 管 
Rj[n]=o2?ay[njas[n] 的 秩 为 1， 但 采样 间隔 过 长 时 ， 对 移动 干扰 源 做 时 间 平 均 产 生 一 个 秩 更 高 的 矩 
阵 ， 这 就 需要 用 更 多 的 特征 矢量 去 张 成 干扰 子 空间 。 只 有 当 采 样 间隔 的 宽度 相对 于 干扰 移动 速率 较 鹤 
时 ， 投 影 算 子 P 才能 有 效 地 投影 这 个 秩 数 较 高 的 干扰 分 量 。 这 个 有 限 窗口 的 估计 带 来 的 影响 通常 叫 
做 “ 子 空间 扩散 ”。 

通过 子 空间 投影 波束 形成 抗 干扰 算法 仿真 得 到 的 在 不 同 输入 干 噪 比 情况 下 单 移 动 干扰 消除 的 干扰 
抑制 比 (Interference Rejection Ratio, IRR) 如 图 1。 

从 图 1 TAER, EPA ERE REAR, SE HA PR EG ri Oe Ee TS lL, 
因为 此 时 R, ,的 主 特征 向 量 更 接近 真实 的 干扰 信号 。 同 时 ， 干 扰 抑 制 比 随 着 样本 数 长 度 的 增加 而 增 
加 ， 直 到 子 空间 扩散 开始 影响 干扰 消除 性 能 ， 这 会 产生 一 个 干扰 抑制 比 的 最 大 值 。 随 着 样本 数 的 继续 
增加 ,干扰 消除 性 能 开始 急剧 恶化 ， 也 就 是 说 干扰 源 的 移动 限制 了 可 得 到 的 干扰 抑制 比 的 最 大 值 。 但 
是 在 较 低 输入 干 噪 比 的 情况 下 ， 由 于 无 法 得 到 干扰 子 空间 的 准确 估计 ， 因 此 无 法 对 干扰 进行 很 好 的 消 
除 ， 所 以 干扰 抑制 比 并 没有 明显 的 变化 。 

通过 子 空间 投影 波束 形成 抗 干扰 算法 仿真 得 到 在 不 同 输入 干 噪 比 情况 下 单 移动 干扰 消除 残余 干扰 
功率 如 图 2。 

从 图 2 可 以 看 到 ， 对 于 一 个 给 定 长 度 的 样本 数 ， 一 个 更 高 的 输入 干 噪 比 带 来 更 低 的 残余 干扰 ， 因 
为 此 时 及 ., 的 主 特征 向 量 更 接近 真实 的 干扰 信号 。 同 时 ， 残 余 干扰 功率 随 着 样本 数 长 度 的 增加 而 减 
小 ， 直 到 子 空间 扩散 开始 影响 干扰 消除 性 能 ， 这 会 产生 一 个 残余 干扰 功率 的 最 小 值 。 随 着 样本 数 的 增 
加 ,干扰 消除 性 能 开始 急剧 恶化 ， 也 就 是 说 干扰 源 的 移动 限制 了 可 得 到 的 残余 干扰 功率 的 最 小 值 。 但 
是 在 较 低 输入 干 噪 比 的 情况 下 ， 由 于 无 法 得 到 干扰 子 空间 的 准确 估计 ， 因 此 无 法 对 干扰 进行 很 好 的 消 


= 除 。 不 过 由 于 输入 的 干 噪 比 很 低 ， 所 以 此 时 的 残余 干扰 功率 也 很 低 。 
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图 1 移动 干扰 源 情 况 下 不 同样 本 数 和 输入 图 2 移动 干扰 源 情况 下 不 同样 本 数 和 
干 嗓 比 下 的 干扰 抑制 比 干 骂 比 下 的 残余 干扰 功率 
Fig.1 IRR dependence on STI length and INR Fig.2 Residual interference power dependence on STI 


level for a moving interferer length and INR level for a moving interferer 
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3 ”基于 多 项 式 模型 的 子 空 间 投影 


由 上 述 分 析 可 知 ， 传 统 的 自 适应 波束 形成 抗 干扰 方法 在 样本 数 有 限 的 情况 下 ， 由 于 协 方差 矩阵 估 
计 存 在 误差 ， 导 致 干 扰 信号 子 空间 估计 精度 受到 限制 ， 从 而 导致 干扰 消除 性 能 下 降 。 另 一 方面 ， 如 果 
增加 样本 数 ， 由 于 干扰 源 运动 导致 干扰 信号 的 子 空 间 扩散 ， 严 重 影响 算法 的 干扰 消除 性 能 。 在 低 干 品 
比 情况 下 ， 由 于 干扰 信号 子 空间 不 准确 ， 传 统 方法 的 干扰 抑制 性 能 不 能 满足 实际 工作 的 需求 。 为 了 解 
决 上 述 问题 ， 本 文通 过 将 时 变 的 R,[n] 与 多 项 式 模型 合并 ， 提 出 了 一 种 子 空间 投影 波束 形成 算法 的 扩 
展 。 这 种 新 的 基于 多 项 式 模 型 的 子 空间 投影 算法 (Polynomial-based Subspace Projection, ，PSP ) 提 高 了 干 
扰 子 空间 的 跟踪 性 能 。 
由 于 干扰 源 的 移动 和 天 线 方向 图 结构 导致 从 R ,到 R,,, 的 变化 是 连续 的 ， 因 此 一 个 适当 阶 数 的 
多 项 式 拟 合 可 行 。 假 设 忽略 射电 天 文 信号 的 影响 ， 这 个 假设 符合 常见 的 射电 天 文 的 观测 场景 ， 因 为 射 
电 天 文 信号 的 强度 过 于 微弱 ， 以 至 需要 长 时 间 的 积分 才 可 以 检测 到 。 
时 间 连 续 的 协 方差 矩阵 尺 (i) ， 给 出 了 它 的 采样 形式 REn]。 引 入 R(t，C) ~R, (CD 作为 时 变 的 干 
> 扰 分 量 R.(1) 的 模型 ， 并 通过 多 项 式 系数 矩阵 C 使 其 参数 化 。 单 干扰 信号 R,y(t)= 02,a,(1) a" (1) IE 
为 1， 因 此 可 以 用 + 阶 的 多 项 式 矢量 构建 时 间 连 续 的 oa lt) 。 


t=[1, t,e, 81"; (7) 
R(t, C)=a(t, C)a"(t, C) = Cu", 
~N 多 项 式 系数 C 的 估计 值 可 以 通过 每 一 个 单独 短 时 积分 窗 的 采样 协 方差 矩阵 的 平方 误差 最 小 准则 得 到 。 


A K-1 NI As 
Ai C =argmin ð, || R, -Rit ©) lI; 
rio C ££=0 (8) 


其 中 , t, HE k MERRIA FEA O ESP PREP ED EEE, (8) ARE TXF O<ASK-1, 
总 有 -Tt 三 7， 在 多 项 式 拟 合 窗 口中 生成 了 个 短 时 积分 窗 。 多 项 式 估 计 的 时 间 尺 度 标准 化 是 任意 
的 ， 因 为 了 不 代表 任何 实际 的 时 间或 以 秒 或 采样 数 为 刻度 的 窗口 宽度 。 在 估计 (6) 式 时 了 作为 一 个 调 
节 参 数 以 限制 时 间 范 围 ， 使 得 高 阶 的 多 项 式 模型 不 会 主导 曲线 走向 ， 并 且 在 这 时 多 个 有 界 的 变化 值 可 
以 用 拟 合 窗口 表示 。 经 过 多 次 仿真 分 析 发 现 ， 将 了 设置 为 0.75 时 得 到 较 好 的 拟 合 性 能 。 

利用 一 致 性 14 12=zrft424 将 (7) 式 代入 (8) 式 ， 并 且 扩 展 到 矩阵 形式 得 到 : 


全 (9) 
f, Si ROR, + CB,C"CB,C"- R,CB,C" - CB,C" R"), 
其 中 ，B =tio BFA (8) SURI IEE Br ASB EE e A RAEE E, A HE A 
法 改善 收敛 时 间 。 
有 一 系列 广泛 的 优化 代码 可 用 来 解决 (9) 式 的 问题 。 第 4 节 的 仿真 结果 通过 MATLAB 的 函数 
fminunc.m 得 到 ， 它 是 一 种 基于 内 部 映射 牛顿 法 的 置信 域 方法 。(9) 式 的 梯度 和 全 局 最 小 解 用 来 确定 
置信 域 问题 并 把 这 个 区 域 限制 在 二 维 子 空间 。 
梯度 : 虽然 多 项 式 的 系数 是 复 值 ， 但 在 平常 的 优化 代码 中 ,它们 必须 作为 独立 的 实 部 和 虚 部 计 
算 。 在 这 种 情况 下 ， 梯 度 需要 改写 为 
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a RefC} 
fara 
VEO) eee 


vec{Re{(CB,C" - R,) CB,}} 


ve( Df) = ‘> | 


1| vee{Im{(CB,C® - R,) CB,} } 
全 局 最 小 解 : 用 EE, 表 示 一 个 除了 (s, 1) 位 置 为 1 以 外 全 0 A EE, ee MR WN 


K K | vec{Re{M,} } vec{- Im{M,} } 
Hel dfs) ~ 1È pea } vec{Re{M,} } | 


: (11) 
M, = E,B,C"CB, + CB,E"CB, + CB,C"E,B, - R,E,B,, 


M, = E,B,C"CB, - CB,E"CB, + t CB,C"E d, - R,E,B, . 


DEI, (ORERE A AMR, RT AE E S, ECAR AE p, A 
(era(t,, C)) (Pal, O) al, Oal, Co i ti 上 ， 而 其 它 的 


we tis = 0 解决 ， 此 
时 有 和 1a(0, C)|,=0。 文 [10] 给 出 的 一 种 参数 约束 的 闪 代 方案 是 使 台 ” 人 1C1,, = 0， 也 就 是 说 


ee N 但 是 这 两 种 方法 
在 系统 中 产生 的 误差 方差 也 一 样 ， 这 些 值 仅仅 依赖 于 相位 的 不 同 。 这 正 是 Wi，C ) 是 如 何 代 入 (7) R 
中 以 得 到 R, C ) 的 ， 因 此 由 于 它 在 优化 代码 中 更 易于 实现 ， 选 择 和 1C1,,=0。 

理论 上 ， 子 空间 投影 算法 可 以 扩展 到 处 理 多 个 干扰 源 的 情况 。 此 时 ，%,(1，C,, ) 是 m 个 干扰 信号 
的 时 变 阵 列 响 应 。(7) 式 改写 为 

RE GE > ct (12) 

在 数值 优化 中 将 梯度 和 全 局 最 小 解 扩展 到 相应 的 C， 中 是 一 项 繁重 却 很 吉 接 的 工作 。 最 大 的 困难 

在 于 为 每 一 个 C， 找 到 一 个 好 的 初始 值 。 因 此 ， 后 面 的 讨论 限定 干扰 维度 为 1 


4 仿真 和 分 析 


在 仿真 中 尽量 保证 仿真 环境 的 设置 与 物理 实 
验 平台 相符 合 。 详 细 的 数值 模型 是 利用 一 个 10 期 望 信号 
阵 元 双 极 化 线性 半 波 长 宽带 的 加 厚 偶 极 子 的 相 控 ty 
阵 馈 源 放置 在 地 平面 上 形成 的 ， 其 结构 如 图 3。 


可 用 的 天 线 带宽 为 300 MHz， 中 心 频 点 为 1 600 TARS 
MHz， 同 时 这 也 作为 仿真 生成 的 窄带 信号 的 中 心 和 NNW 


频率 。 数 据 的 采样 速率 为 1. 25 Msamp/s， 带 宽 
为 450 kHz, 
以 调制 C/A 码 的 北斗 导航 信号 为 干扰 信号 。 辅助 天 线 YN 
短 时 积分 窗 序 列 中 的 每 一 个 窗口 都 通过 2 ~27 人 
个 采样 点 计算 ， 并 且 所 有 的 短 时 积分 窗 使 用 采样 ae 
Fig.3 The Schematic diagram of simulation 


估计 校准 的 噪声 协 方差 矩阵 尺 进行 预 白化 。 同 antenna structs 
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时 每 一 个 相位 正确 的 特征 矢量 进行 多 项 式 拟 合 。 利 用 (5) 式 得 到 的 静态 (无 干扰 ) 波 柬 形成 权 矢 量 
w= Roa, 是 最 大 灵敏 度 的 波束 ， 此 时 有 用 信号 保持 在 天 线 的 视 轴 上 。 

图 4 给 出 了 在 每 一 个 短 时 积分 窗 的 协 方差 矩阵 中 使 用 多 项 式 模型 拟 合 与 传统 的 子 空间 投影 算法 的 
性 能 分 析 对 比 。 从 图 中 可 以 看 出 ， 子 空间 投影 算法 对 短 时 积分 的 窗口 宽度 有 更 高 的 敏感 性 ， 但 是 子 空 
间 投 影 算法 在 整个 短 时 积分 窗 的 采样 窗口 的 性 能 优 于 子 空间 投影 算法 。 在 很 短 的 短 时 积分 窗 采样 窗口 
情况 下 ， 子 空间 投影 算法 的 性 能 急剧 恶化 ， 这 是 因为 此 时 的 采样 协 方差 矩阵 有 太 多 的 估计 误差 ， 不 能 
得 到 一 个 良好 的 特征 向 量 。 而 子 空间 投影 算法 在 较 短 的 短 时 积分 窗口 宽度 的 情况 下 可 以 得 到 与 较 长 的 
短 时 积分 窗口 宽度 近似 的 干扰 消除 性 能 。 图 3 和 图 1 的 仿真 环境 都 是 在 相同 干 噪 比 和 3° /s 的 干扰 移 
动 速率 下 进行 。 在 图 3 中 ，PSP-8 是 指 特征 矢量 是 用 一 个 8 阶 多 项 式 拟 合 的 ，PSP-12 是 指 特征 矢量 是 
用 一 个 12 阶 多 项 式 拟 合 的 。 


-20 lutput INR -12dB, 3 /sec motion -24 lutput INR 0dB, 3°/sec motion 
25 —— SP -28 一 6 一 SP 
s Š —e— PSP-8 
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图 4 不 同 的 干 噪 比 和 干扰 移动 速率 下 的 子 空间 投影 算法 和 基于 多 项 式 的 子 空间 投影 算法 的 性 能 。 
(a) INR=-12 dB; (b) INR=0 dB; (c) INR= 12 dB 
Fig.4 Performance of SP and PSP for various input INRs and motion rates 


5 结 论 


本 文 的 研究 要 点 在 于 : (1) 详 细 分 析 了 子 空间 投影 算法 的 性 能 特点 ， 并 研究 了 限制 干扰 消除 性 能 
的 因素 ， 包 括 采 样 估计 误差 和 子 空间 扩散 ;，(2) 提 出 了 一 种 多 项 式 拟 合 的 方法 以 更 好 地 构造 由 干扰 移 
动 引起 的 时 变 特性 。 虽 然 本 文 是 使 用 射电 天 文 的 相 控 阵 僻 源 仿 真 和 分 析 算 法 的 有 效 性 ， 但 对 于 其 它 的 
移动 干扰 环境 下 的 自 适 应 阵列 干扰 消除 有 广泛 的 适用 性 。 

更 进一步 的 工作 包括 把 算法 应 用 到 多 干扰 消除 和 处 理由 多 径 散 射 带 来 的 秩 的 增加 。 本 文 的 实验 数 
据 是 以 秒 为 单位 进行 观测 的 结果 ， 但 是 实际 的 射电 天 文 观 测 中 可 能 需要 以 小 时 为 单位 进行 观测 。 因 
此 ， 实 时 处 理 需要 考虑 该 方 法 是 否 可 以 在 实际 的 天 文 观 测 中 应 用 。 实 时 处 理 的 问题 包括 快速 多 项 式 参 
数 估 计 和 波束 形成 中 的 自动 参数 估计 问题 ， 因 此 需要 根据 不 同 的 数据 环境 选择 不 同 的 多 项 式 阶 数 。 
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Abstract; Radio astronomical data are increasingly corrupted by human telecommunication activities. 
Therefore, Radio Frequency Interference ( RFI) mitigation becomes an important step in the data processing 
flow. This paper considers the problem of adaptive array processing for interference canceling to drive nulls in 
a moving interference environment. As to the telescopes based on antenna array, RFI mitigation can be 
conducted in spatial domain using the sampled covariance matrix. In many practical scenarios, the achievable 
null depth is limited by covariance matrix estimation error which leads to poor identification of the interference 
subspace. We pay attention to the particularly troublesome cases of relatively rapid interference motion. A 
polynomial-based model is incorporated in the proposed algorithm to track changes in the array covariance 
matrix over time, mitigate interference subspace estimation errors, and improve canceler performance. 
Performance for conventional subspace projection (SP) is compared with polynomial-augmented SP using 
simulation. The simulation results validate the effectiveness of the proposed algorithm. 
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